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'-H ! Abstract 

Oh, Lattice proposals for a nonperturbative formulation of the Standard Model easily 

^ ■ lead to a global U(l) symmetry corresponding to exactly conserved fermion number. 

The absence of an anomaly in the fermion current would then appear to inhibit 
anomalous processes, such as electroweak baryogenesis in the early universe. One way 
^ ! to circumvent this problem is to formulate the theory such that this U(l) symmetry 

■ is explicitly broken. However we argue that in the framework of spectral flow, fermion 

creation and annihilation still in fact occurs, despite the exact fermion number 
conservation. The crucial observation is that fermions are excitations relative to the 
vacuum, at the surface of the Dirac sea. The exact global U(l) symmetry prohibits 
a state from changing its fermion number during time evolution, however nothing 
prevents the fermionic ground state from doing so. We illustrate our reasoning with 
a model in two dimensions which has axial-vector couplings, first using a sharp 
momentum cutoff, then using the lattice regulator with staggered fermions. The 
difference in fermion number between the time evolved state and the ground state is 
indeed in agreement with the anomaly. A study of the vacuum energy shows that the 
perturbative counterterm needed for restoration of gauge invariance is insufficient 
in a nonperturbative setting. For reference we also study a closely related model 
with vector couplings, the Schwinger model, and we examine the emergence of the 
^-vacuum structure of both theories. 
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1 Introduction 

A tractable nonperturbative formalism for the analysis of cliiral gauge theories is 
an outstanding problem in modern field theory. Among the difficulties it has been 
argued that lattice formulations (for reviews see 0) or indeed any regularization, 
cannot yield the correct physics if it implies an exact global symmetry which should 
be broken by anomalies 0, ^. For example, such models would be unable to 
describe anomalous fermion number violation in the Standard Model, which for 
example, plays an important role in current explanations of the observed baryon 
asymmetry in the universe. 

A lattice regularized version of the Standard Model is of course supposed to 
describe these effects correctly. However typical forms of the lattice fermion ac- 
tion, Sf = —J^xyi^x^xyi^y havc an exact global U(l) invariance ip ~^ exp(itt; )■?/', 
ip — >■ exp{—iuj)ilj and since the lattice fermion measure is generally also invari- 
ant, one expects this global symmetry to be accompanied by exact fermion number 
conservation, suggesting that such models cannot display the anomalous fermion 
production of the Standard Model. We shall call this the electroweak U(l) problem. 

Some recent discussions have focused on the question of heavy fermion decou- 
pling ll], ^, 1^, 1^, but the U(l) problem is more general. We also remark that its 
resolution should not depend on alternative definitions of the fermion number cur- 
rent, constructed such that their divergence has the right anomaly (see e.g. @, 0). 
Once the action and measure in the path integral are defined nonperturbatively, the 
stage is set, and we should just have to see what physics emerges. 

One way out would be to construct lattice actions such that they have no un- 
wanted exact global symmetries ^, ||, ^ |T0[ and it is reasonable to assume that 
the violation of the global symmetry turns in the scaling region into the desired 
anomaly [|TT| . 



Here we want to reconsider this standard lore and study the effects of the un- 
wanted global U(l) symmetry in general terms. We examine the spectral fiow to see 
whether creation or annihilation of fermions, i.e. fermionic particles, takes place in 
a simple model in 1+1 dimensions which has the exact global U(l) symmetry corre- 
sponding to fermion number conservation. We study the spectrum of the fermionic 
Hamiltonian in an external gauge field which changes slowly in time, the idea being 
that we consider time slices of an instanton-like configuration, which starts as a vac- 
uum gauge field, goes through a sphaleron configuration and ends up as a different 
vacuum gauge field related to the first by a topologically non-trivial ('large') gauge 
transformation. If the change with time is very slow, we can make use of the adi- 



abatic theorem [|12| and deduce the evolution of the states by continuity, following 
the eigenstates of the time dependent Hamiltonian. We can then check explicitly 
whether an initial fermionic vacuum state has evolved into a state with particles 
produced by the gauge field in accordance with the anomaly. Particle annihilation 
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proceeds of course by the inverse process. For a recent discussion of anomalies and 
spectral flow see ref. [|l3l; an earlier spectral flow analysis in lattice gauge theory 



was given in ref. [M and more recently in the context of domain wall fermions [IT 



in ref. [16 



Based on the spectral flow analysis we argue that fermion production occurs 
despite the exact fermion number conservation. The crucial observation is that 
fermions are excitations relative to the vacuum, at the surface of the Dirac sea, which 
are defined as the ground states in external fields that are 'pure gauge'. The exact 
global U(l) symmetry prohibits a state from changing its fermion number during 
time evolution, however nothing prevents the fermionic ground state, the state of 
lowest energy which depends on the adiabatic external gauge field, from changing 
its fermion number. We show that the difference in fermion number between the 
time evolved state and the ground state after the sphaleron transition is indeed in 
agreement with the anomaly. 

Our arguments do not really depend on lattice regularization. Still, we would 
like to indicate how they might fit into a possible lattice formulation of chiral gauge 
theories, taking into account the experience gained from investigations over the 
past years showing the failure of many promising proposals [Q. Recall an often 
followed practice in the continuum: the fermions are integrated out first and from 
the resulting determinant a well defined effective action is constructed for the bosonic 
fields. This can be done in a variety of ways; in general it requires the addition of 
counterterms to get a finite answer upon removing the fermionic regulator, assuming 
the fermions are in an anomaly free representation of the gauge group. Subsequently, 
the integration over the bosonic fields is to be performed. One first removes the 
fermion regulator and then the boson regulator. We see no reason in principle why 
such a practise cannot be followed on the lattice. It means dealing with lattice 
fermions in smooth external gauge fields, a relatively simple situation for which any 
lattice fermion method should work. Counterterms are usually needed to restore 
(chiral) gauge invariance, and with a gauge invariant effective action the introduction 
of ghost fields Jl^ may be avoided as in non-chiral lattice gauge theory (the problem 
noted in |T^ need not apply). Taking the fermion lattice distance all the way to 



zero first ('the desperate's method' 0, |T9|), may not be necessary if the violation 



of gauge invariance can be controlled so that gauge symmetry restoration may be 
invoked |^ . Obviously, much work still remains to be done in this direction, and our 
present study of the subtleties of anomalous fermion production seems a necessary 
step, which is also of interest in its own right. A preliminary account of this work 
has already been presented in ref. |^ |. 

The outline of the paper is as follows: In Sect. 2 we introduce the two-dimensional 
massless models, the Schwinger model and its equivalent axial version (axial QED2), 
which we shall use in this paper to illustrate our reasoning. In Sect. 3 we first discuss 
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the spectral flow of these models with a sharp momentum cutoff on the number of 
modes. We argue that fermion creation and annihilation in axial QED2 is possible 
in spite of the exact U(l) symmetry. The sharp momentum cutoff is unsatisfactory 
because it is both non-local and gauge variant. Local lattice versions of the two 
models are introduced in Sect. 4, using staggered fermions. The lattice axial QED2 
model lacks gauge invariance, which can be restored for smooth external gauge fields 
with a mass counterterm for the gauge field. We then study the spectral fiow in 
both models to see how the properties of the well understood vector model compare 
with those of the axial model. Fermion creation and annihilation takes place in 
the same way as in the toy models with the sharp momentum cutoff. There are 
furthermore interesting aspects to the way the counterterm restores gauge invariance 
to the ground state energy in the axial model, and to the way the Dirac sea acquires 
a gauge invariant bottom in the vector model. We summarize our conclusions in 
sect. 5. For clarity of presentation we have delegated many details of the staggered 
fermion formalism and the calculation of the energy spectrum to appendices A-C. 

2 Axial and vector QED2 in the continuum 

The massless axial QED2 model is given in the continuum by the following (real 
time) action, 

S ^ - J d'x^F,,F^^ + Sf, (2.1) 

Sf = - J d^xi^Yidi. + tA^'y,)^;, (2.2) 

where 7''^ = 71 = (Ti, 7° = —70 = — i(J2, 75 = Cs- We take to be an external 
gauge field; only the fermion fields are quantized. The model is invariant under local 
axial gauge transformations: ipi^) ~^ exp(ia;(a;)75)-?/'(a;), ipi^x) ipi^x) exp^iuj^x)^^) , 
Ai^{x) — >• Ai^{x) + d^u!{x), with gauge current = iipj'^j^ip. The action is fur- 
thermore invariant under the global U(l) symmetry ip{x) —>■ exp{ia)ip{x), iplx) — > 
ip{x) exp(— ia), with a corresponding vector current — iip^^ip. According to stan- 
dard lore the gauge current, upon quantization, has to be divergence free, 9^j5 = 0, 
while the divergence of the vector current becomes anomalous 

d,f = -2g , q= ^e^'F^. ^ d,C>^ , (2.3) 

where C^{x) — {l/27r)€^''Aiy{x) is the Chern-Simons current and q{x) is the topo- 
logical charge density (eoi = +1); the topological charge is defined by / (Px q{x). 
Fermion and Chern-Simons numbers are given by the spatial integrals over the time 
components of the corresponding currents, Q = J dx^ j'^ix), C = J dx^ C^{x). For 
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later purposes we also introduce here the axial charge Q5 = J dx^ jiix) which is 
exactly conserved in this model. From eq. (f^.3D one can deduce the equation 



Q{t) - Q(0) = -2{C{t) - C(0)) = -2 dt J dxq , x = x\ t = x\ (2.4) 

which relates a change in the Chern-Simons number to a change in fermion number. 
The aim of this paper is to investigate whether, in those models which have the 
exact global U(l) symmetry and where Q is thus conserved, a change in C due to a 
sphaleron transition can still give rise to a change in the number of fermions in the 
spirit of eq. (|2.4|). 

By charge conjugation of the right-handed fermion fields, 

^« = (?^Cf, = -(CV^i)^ , ^L = V^1, ?L = ^L (2.5) 

with C the charge conjugation matrix, or alternatively to only the left-handed fields, 
one can show that axial QED2 is equivalent to the massless Schwinger model (vector 
QED2), 

S'f = -I d^x'^'rid^ - iA^)ij' . (2.6) 

For the moment we will label the quantities in the vector model with a prime for 
clarity. The local axial gauge symmetry of the action ( |2.2|) has turned, after the 
transformation ( p.5|) , into a local vector gauge symmetry for the ip' fields, il)'{x) — > 
e'X]i{iuj{x))ip' {x) , ip (x) ip (x) exp{—iLj{x)), and vice versa the global vector invari- 
ance into a global axial invariance ip'{x) exp{ia'j^)il!' (x) , ip (x) ^ ip (x) exp(ia75). 
The corresponding currents and charges transform under (|2.5|) into each other. 



Q = -Q',, Qs = -Q' . (2.7) 



We notice that a mass or Yukawa term added to the actions (|2.2| ) and ( p.6|) would 
transform under ( |2.5| ) into a corresponding Majorana mass or Yukawa term. In this 
paper we will not include these terms. For a discussion of the spectral flow in QED 
in presence of a bare mass term see ref. |n . 



In the following we use a finite spatial extent of length L, with periodic boundary 
conditions for the gauge field and antiperiodic boundary conditions for the fermion 
fields. We consider gauge fields with vanishing time component, Ao = 0, and choose 
the Coulomb gauge specified by Ai{x,t) = A{t)] then the time dependence is ex- 
pressed by A-dependence. In this gauge the Chern-Simons number C = —AL/2n, 
and going through an instanton-like configuration means that C changes by one 
unit. Values A = 2TTk/L with integer k are 'pure gauge', since for these values we 
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can write A = Qidx^*, with Q = exp{i27Tkx/L) a periodic gauge transformation 
with winding number k. When A = 2TTk/L the Chern-Simons number takes integer 
values C = —k. An example of such gauge fields is given by the configuration with 
constant electric field / = Fqi, 

Ao = 0, A^{x,t) = A{t) = -ft. (2.8) 

For / = 27i/LT we have C(0) = 0, C{T) = 1, and the adiabatic limit corresponds 
to T — > oo. 

The ground state \0,A) of the fermionic Hamiltonian H{A) is by definition the 
state with lowest energy. When A is pure gauge, these ground states are to be iden- 
tified with the vacuum. In Hilbert space, gauge transformations on A induce unitary 
transformations, and we have to identify states related by gauge transformations, 
hence also all vacua with differing integer Chern-Simons numbers. 

A quantity which we will use in our analysis is the ground state energy Eq{A), 
H[A)\0, A) = Eo{A)\0, A). It can be obtained from the partition function at inverse 
temperature f3 in the limit (3 ^ oo 

Eo{A) = -^InTre-'^^^^), (3 ^ oo , (2.9) 

or in terms of the euclidean effective action in a space-time volume L x /3, 

S'eff(A) = In/ Di^Di: exp(^7^'), (2.10) 
^ -f3Eo{A), (3^oo. (2.11) 

We shall use the euclidean formalism for our lattice models. The quantity Eq{A) 
can be evaluated as (see for example |^ and references therein) 

Eq(A) = — y ^— ^ cos(nLA) + const. (2.12) 

The sum is clearly periodic in A with period 2it/L and evaluates, in each period, to 
a quadratic potential. 

Eo{A)-EoiO) = ^A'', -l<A<+y, mod(27r/L). (2.13) 



From ( p.l3|) we can read off the dynamically generated mass of the Schwinger Model, 
= e^/vr (the factor is due to the term in (|2.1| )). For a treatment of the 
Schwinger model on a circle in the hamiltonian formalism, see refs. p3 |. 
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3 A model with a sharp momentum cutoff 

Before giving a full treatment using lattice fermions we first consider a model of 
fermions regulated by a sharp cut-off in momenta A. Such a regulator captures 
essential features we wish to study in lattice regularization, namely, exact global 
invariance and restriction to a finite number of modes for the ip field. It is not gauge 
invariant, which it shares with the class of chiral lattice models we have in mind, and 
moreover it is non-local, whereas the lattice models are local. The locality of the 
lattice models allows for restoration of gauge invariance by a local mass counterterm 
oc / (PxAf.A^'. 

From the above actions (|2.2| ) and (p.6|) , the axial QED2 hamiltonian H and the 
equivalent vector QED2 hamiltonian H' are given by 

H = J T.[^Up)Mp)iP + i^lip)Mp)i-P + A)] - AN/2 , (3.1) 

H' = jT.[^Up)^'r{p){p-A)+^2{p)^M{-P + A)], (3.2) 
^ P 

where the summation is over momenta p = [n — l/2)27r/L G [—A, A] (antiperiodic 
boundary conditions for the fermions), is the total number of modes, and L 
and R denote the left-and right-handed projections 75 — >■ =f1. The Fourier modes 
ip{p) = Jq dx exp{—ipx)ilj{x) have the commutation relations {ip{p),i^^{q)} = LSp^q. 
In axial QED2 the spatial component of the gauge field couples to fermion number 

Q = J dx^[ip\x),'ilj{x)] 

= j&Up)Mp)+i^lip)Mp)]-N/2. (3.3) 
^ p 

Note that upon quantization we have anti-symmetrized the field operators in Q, 
which is necessary for it to change sign under charge conjugation. Charge conjuga- 
tion ( |2.5D on the right handed fields only, relating the axial and vector model, is given 
by i^R = i^R, 4^R = 4^R- Denoting this conjugation by Cr, we have H' = CrHCI,, so 
that H' has the same spectrum as H. 

In these models regularized by a sharp momentum cutoff, there are two conserved 
charges: both Q and 

Q5 = / dx^\xh,^{x) = J Y.['^Up)Mp) - ^Lip)Mp)] (3.4) 



commute with H. In the vector version the two charges are given by Q' = —Q5, 
Q'^ = -Q (cf. eq. ^). 
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The spectral flow of the eigenmode energies e in the axial model is very simple. 
Disregarding the mode-independent term —AN/2 we have 



(3.5) 



where the chirality x — is the eigenvalue of 75 (x = +1 (—1) for the R (L) 
modes). Each mode's energy increases linearly with A. In the vector version 

e'^{p,A) = x{p-A) , (3.6) 

and the it! (L) modes decrease (increase) hnearly with A as shown in fig. 1. 



Axial 



Vector 




0.2 0.4 0.6 0.8 
hA/2n 




0.2 0.4 0.6 
LA/27r 



Figure 1: Occupation of the ground state \0,A) (solid and dashed) and the adiabati- 
cally time evolved state A) (solid and dots) in the momentum cutoff models. 



Consider now the flow of the state which starts out as the vacuum state 

at A = 0, i.e. the Dirac sea with all negative energy modes occupied. Each mode is 
doubly degenerate (L and R). In the axial model the flow upwards with A and at 
A = 271 /L they have taken the place of their predecessor, except for the two modes 
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starting out as e = —tt/L and ending up at e = +11/ L. The quantum numbers Q 
and Q5 are conserved, both are zero for all states A), as for the initial vacuum 
state |0,0). The energy Eq, of the state l^',^) also happens to be constant in this 
regularization: each mode shifts upwards by A and the N/2 modes contribute AN/2, 
which is compensated by the mode-independent term —AN/2. However, the ground 
state |0, A) looses two occupied modes half way at A = n/L, when e = —tt/L + A 
crosses zero. The final vacuum energy Eq at A = 27t/L differs from the initial one, 
due to the lack of gauge invariance of the sharp momentum cutoff. 

The important point is that the final ground state has two occupied states (L 
and R) less than the initial ground state, so the final vacuum quantum numbers are 
Q = —2 and Q5 = 0. Hence, 

AQ = {Q,^ — (5o)final — (Q-I- — Qo)mitial = +2 , (3.7) 

just as expected from the anomaly in the vector current, AQ = —2AC. In the 
vector model the L modes move upwards and the R modes move downwards, with 
(Q5,*)initiai = (Q5,*)finai = 0. Midway, the ground state looses an L mode and gains 
an R mode, such that (Q5,o)initiai = 0, (Q5,o)finai = +2 and AQ' = 0, AQ'^ = -2 = 
2 AC. In the axial model we have the creation of two particles (L and R), whereas 
in the vector model we have the creation of a particle (L) and an antiparticle (R) 
(an i?-hole). 

In both models the final energy difference of l^',^) with [0,^4) is two units 
tt/L, which corresponds to the creation of the two particles (using the generic term 
'particle' also for the antiparticles). It is amusing that this can be interpreted in 
classical terms, as the work which is provided by the external field in creating two 
fermions [|l^]. The time component of the Lorentz force dp^/dt = -^F^^^dx^ /dt on a 



particle of charge ±1 can be integrated for the gauge field of constant field strength 
( p.8|) . Using dx^/dt = ±1 for the massless particles, and taking into account that 
they appear midway at t = T/2, this gives a final energy per particle p^{T) = 
dp^/ dt = it/ L. 

Obviously, the ground state energy in these cutoff models is very different from 
the continuum form ( ^.131 ). It is given by Eo{A) — Eq{0) = for < A < ti/ L and 
Eo{A) - Eo{0) = -2{A- -k/L) for vr/L < A < 27: /L. A 'counterterm' piece- wise 
linear in A would be required to rectify the situation. This is at variance with the 
expectation that a simple quadratic mass counterterm oc / d'^x A^A^ should be suffi- 
cient to restore gauge invariance. The sharp momentum cutoff is also unsatisfactory 
because it is non-local. The local lattice models to which we shall turn in the next 
section are well behaved in this respect. 
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4 Level shifting in the lattice models with stag- 
gered fermions 

Our lattice axial and vector QED2 models are formulated in terms of staggered 
fermion fields. In order to keep the technicalities to a minimum, we have collected the 
details of this formulation to the appendices and here we just record the essentials. 
We use the Euclidean lattice formulation with the partition functions given by 

^ = ([[dXxdXx) expSf, (4.1) 

X 

with 

Sf = -Yv^.xl{XxUf,xXx+^ - Xx+f,U*xXx) , (4.2) 



Xfl 



for vector QED2 and 



Sf = -Y.\Y.'^^^^^^^V^^MXxXx+f,-Xx+^,Xx) 

X { fl 

+ ^ sin A^a; e^jyT/jyj;. 2 (XxXx+s + Xx+pXx) r + '^o^/ix 

flU ) Xfl 



(4.3) 



for axial QED2. We use lattice units, in which the lattice distance a = 1. U^^x = 
exp(— iA^a:) is the lattice gauge field, the rj^x are the usual staggered fermion sign 
factors representing the gamma matrices, r/i^ = 1, r]2x = (~1)^^ (72 = ^7°, ^12 = +!)• 
The last term in ( |4.3|) is a counterterm. The vector model is gauge invariant; the 
axial model is not, however invariance is restored in the scaling region by the addition 
of this counterterm. The derivation of the continuum models ( |2.6| ) and ( p.2|) from 
the staggered fermion actions (^4.2|) and ( [4.3|) is sketched in appendix A, where 
we give also the perturbative result for the coefficient r of the counterterm. A 
nonperturbative extension of the counterterm is discussed below. 

The exact U(l) phase invariance of the actions and measure in the path integrals 
above is obvious. For the vector model this is just the global limit of its local gauge 
invariance. For the axial model we have an exact global U{1) symmetry leading 
to fermion number conservation. To derive the divergence relation of the current 
associated with this symmetry we replace Xx Q^V>ii^x)Xx and Xx ~^ 6xp(— zci;^.)x^ 



in the path integral ( |4.1| ) and note that it is invariant under this transformation. 
After collecting terms linear in ujx we find from the fact that their coefficients must 
vanish, the exact current divergence relation on the lattice, 

EOmx-Jm-a)x = 0' (4-4) 
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where the brackets denote the average over the fermion fields. For the vector theory 
the current is given by 

jt,x = iv,^x\{XxUi,xXx+[i + Xx+fiU*^xXx) , (4.5) 
while for the axial theory 

3^lx — ^ '^'^^ ■^^lx'^^lX2i.XxXx+^l ~\~ Xx+jiXx^ 

+ i sin A^^ e^^T]^^ \ iXxXx+fi - Xx+^Xx) ■ (4.6) 
The current in the vector theory ( [4.5|) is of course gauge invariant, but this is lacking 



in the axial theory ( [4.6|) . Presumably the combination j^x+'^iC^x = jiix-\-{i/'^)^iiv^ui 
which has the correct euclidean anomaly, will become gauge invariant in the scaling 
region and can be compared with the continuum current in ( p.3| ). 

The above staggered fermion models actually describe axial or vector QED2 with 
two fiavors. The corresponding SU(2)ixSU(2)/j fiavor symmetry is reduced on the 
lattice to the exact global U(l) symmetry Xx — ^ exp(zcc;ea;)x2:5 Xx ~^ ^w{'^^^x)Xxy 
= (— 1)^^"*"^^, which corresponds to a fiavor non-singlet axial transformation in the 



scaling region ||30[ . The charge corresponding to this symmetry is given in appendix 
B. As long as the gauge field is external, the two fiavors are an inessential extension 
of the previous models. 

Other currents may be constructed which reduce in the scaling region to the 
expected continuum currents. These are somewhat arbitrary because they do not 
correspond to symmetries of the action, although natural choices exist 0. It has 
been shown in ref. that the natural choice for the axial current in the gauge in- 
variant vector QED2, for example, has the expected anomalous divergence equation; 
see for t/(4) x t/(4) currents in staggered fermion QCD. 

The gauge potential A^^. is taken to be independent of euclidean ('imaginary') 
time X2- It depends parametrically on the real time t. From the euclidean path 
integrals above we can derive the transfer operator T in the standard way (cf. ap- 
pendices B and C), Z = TrT^/^, where /3/2 is the number of pairs of time slices 
and f3 the inverse temperature. The hamiltonian H is defined in terms of T by 
T = exp{—2H), such that Z = Ti exp(— /5if). 

For the simple case Ai = A(t) (i.e. independent of xi and X2), A2 = 0, we can 
obtain the spectrum of the transfer operator explicitly, and from this the energy 
levels of the associated hamiltonian as a function of t — the spectral fiow. From 
the calculations outlined in appendices A-C we find the following formulas for the 
energies. For vector QED2: 

ex(P'^) = eox(p-^), 



^Oxip) = ln(x sinp -|- ysin p -|- 1) , or sinh eo^^- = x sin p, (4.7) 

27r vr L L 

P = n---, n = -- + l,...,-. 
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Here the ^oxip) are the free staggered fermion energies and x = =tl is the chirahty 
of the mode, the eigenvalue of 75. Each mode is doubly degenerate because of the 
two flavors. The eigenvalues of the hamiltonian, 

E = L\n2+ ^ ex{p,A), (4.8) 

occupied 
modes 

consist of the sum of the occupied mode energies plus an overall constant Lin 2, 
which could be avoided by redefining the fermion measure in the partition function. 
For axial QED2: 

/ N -, / 1 + sin y4 . , . 

e^{p,A) = eo-x{p)+ln^y-—j , (4.9) 

where we assume cos A > (below we shall furthermore restrict ourselves to \A\ < 
7r/4). In this case 

E = Lln2-Lln(l + sinA)+ ^ e^{p,A), (4.10) 

occupied 
modes 

where E does not include the counterterm. The term —L ln(l+sin A) is the analogue 
of the 'charge term' —AN/2 in (|3.1| ), generalized to two flavors (in lattice units 
L = Na = N). 

We have plotted in fig. ^ the energies ex{p,A) as a function of AL/2tt in the 
interval between and 1, for a lattice of size L = 32. Each energy eigenvalue in 
(|4.7|) and ( [4. 91) appears twice because of the two flavors. For momenta p and gauge 
fields A which are small compared to the cutoff 7r/2, ( [4. 7] ) and ( |4.9| ) reduce to the 
linearly spaced modes in formulas ( |3.6|) and ( p.5| ). The spectra close to the surface 
of the Dirac sea, for small A, are as in the continuum. For the modes near the 
maximal negative or positive energies we clearly see regularization effects. But most 
importantly, the interpretation of anomalous fermion creation and annihilation for 
the spectra in fig. is just as in the previous section: a given quantum state respects 
the global U{1) invariance of the action, while the vacuum adjusts, producing or 
absorbing the charge in concordance with the anomaly. 

There are furthermore some interesting details which are worth mentioning. In 
the vector case the two spectra describing the vacuum at A = and A = 2tt/L are 
identical, in accordance with gauge invariance. As in the sharp cut-off model, the 
L-modes move upwards, whereas the i?-modes move downwards. This however does 
not hold for the upper- and lower-most modes. For example, the maximal negative 
energy mode appears to start out right handed A = (going down) while ending 
up left handed (going up). This effective chirality flip occurs because the staggered 
fermion analogue of 75 (p — A) is 75 sin(j9 — A) the slope of which, {d/dA) 75 sin(j9 — 
A) = —75 cos (p — A), changes sign at A = n/L when the momentum is minimal. 
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Figure 2: Spectral flow for the axial and vector models with staggered fermions. The 
energies e^{p,A) are shown as function of LA/2TT for L = 32. Ground state \0,A) 
and A) are represented as in fig. 1 



p = — 7r/2 + tt/L (at the bottom of the Dirac sea). From a more general point of 
view we expect in coupled staggered fermion theories to be able to recover 75 only 
in the scaling region, i.e. for sufficiently small p. Far from the scaling regime the 
labels 'L' and 'i?' may lose their usual meaning. In the present situation, however, 
the perfectly smooth (constant) A allows for a more detailed identification of 75 also 
outside the scaling region, which leads to the chirality fiip in the minimum energy 
modes at A = n/L. 

The chirality of an eigenmode may be viewed as the eigenvalue of an axial charge 
Qs, as defined in ( |B.27| ), which commutes with the transfer operator. This Q5 is 
the analogue of / dxip^'y^ip — 2C in continuum considerations and it is not gauge 
invariant. On the other hand the natural gauge invariant but not conserved Q5 
constructed from the current j^s = —iSfj^uju gives only sensible chiralities in the 
scaling region (cf. 



The behavior under large gauge transformations is also interesting. Although the 
spectrum is invariant, e.g. comparing A = and A = 27r/L, the eigenmodes shift 
into each other, iPl,r{p) V'-L,-R(p~27r/L), except for the last modes with minimum 
momentum pmin = — '/r/2 + tt/L which turn into the Pmax = 7r/2 — vr/L modes while 
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flipping chirality and flavor. This follows from the transformation properties (B.22) 



of the eigenmodes. The Dirac sea vacua have all positive momentum L-modes and 
negative momentum i?-modes fllled, such that |0); = |0,y4 = 2ttI/L) transforms 
under large gauge transformations fife(x) = exp(z27rx/c/L) as 

f^fe|0)i = \0)k+u (4.11) 

assuming a suitable choice of phases. The action of on these states is given by 

QslO)^ = 2fc|0)fc, (4.12) 

which also expresses the non-gauge invariance of this operator. 

The vector QED2 ground state energy Eq{A) is periodic in A with period 2t^/L, 
which expresses invariance under large gauge transformations. In the continuum 
limit it reduces to the quadratic form ( |2.13|) . The subtle regularization effects in 
the spectrum away from the scaling region concord with gauge invariance and the 
A dependence of Eq is a surface effect of the Dirac sea (cf. eq. ( |A.16| )). In physical 
units, the relative discretization errors in the A dependence are of order a?' A?. In 
the inflnite volume limit the ground state energy density £0 deflned by 

£,{A) = hm (4.13) 

becomes independent of A, £o{A) = £oiO), since the variation in Eo{A) is of order 
1/L (cf. (^)). 

For the axial model flg. § shows that the breaking of gauge invariance has the 
effect that the spectra describing the vacua at A = 27r/L and A = differ from each 
other, especially for high momentum modes which are far away from the surface 
of the Dirac sea. The regularization dependence of the axial spectrum is also a 
subtle effect, which is just right so that the quadratic counterterm can restore gauge 
invariance to the ground state energy in the scaling region. In flg. ^ we have plotted 
the dimensionless function L{Eq{A) — Eq{0)), without the counterterm, with the 
counterterm — rL^A^/2, and with a more reflned subtraction in terms of the energy 
density of the ground state, —L'^£o{A). The latter subtraction is equivalent to the 
perturbative counterterm in the scaling region since £o{A) —£q{0) = tA"^ /2 + 0{A^) 
(cf. appendix A), but leads clearly to better overall behavior for the renormalized 
ground state energy. (The size of the scaling region will be discussed shortly.) 

Fig. ^ shows a close-up. We clearly see how the counterterm restores gauge 
invariance in the scaling region, which for a lattice of 32 sites appears to be given by 
\LA\/2n<l. The cusps at LA/27C = ±1/2 correspond to the energy barrier between 
the vacua Eq{tt / L) — -Eo(O) = / L (for two flavors). In the region —1 < LA/2t{ < 1 
containing the flrst three minima, the vector and axial ground state energies are 
very close. 
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LAEo(A) 




Figure 3: The scaled energy difference LAEq = L{Eo{A) — Eo{0)) for the axial 
model as a function of AL/2tt, on a lattice with N = 32. Shown are the unrenor- 
malized energy: dashed; with perturbative counterterm —tL?'A^/2 added: dotted; 
with L'^{£q{A) — £^o(0)) subtracted: solid. 



Notice that the abscissas in figs. ^ and ^ are hmited to \A\ < 7r/4; the reason 
is that at y4 = ±7r/4 (up to corrections of order 1/L) even the lowest /highest lying 
mode has been pulled out/pushed into the Dirac sea: e^(p, A) = for p = =F7r/2, 
A = ±7r/4. When the Dirac sea evaporates {A ?s vr/4) or the Dirac sky condenses 
{A ^ — 7r/4), the physics clearly ceases to make sense, and we have to restrict the 
allowed values of A in the axial model to \A\ < tt/4. 

We also see from figs. ^ and ^ that the quadratic counterterm is much too 
strong in the large field region, outside the scaling region. With dynamical gauge 
fields it would be bad to have the minima of the classical energy lying outside 
the scaling region and away from A = 0. To bring the minimum of the energy 
into the scaling region we must add suitable higher order terms to the action, e.g. 
— / (i^xr'a^(A^A^)^, with r' > 0. In the static case the best choice is perhaps to 
subtract L£o{A) from the energy. 

We have checked numerically using a variety of lattice sizes that for given LA, the 
relative discretization errors are again O(a^A^). However, in contrast to the vector 
case, convergence is not periodic in LA. We can estimate the size of the scaling region 
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Figure 4: As in fig. on a smaller scale (clashed; unrenormalized energy; dot- 
ted; with perturbative counterterm —tL? A? jl added; solid; vjith L^(£^o(^) ~ ^o(O)) 
subtracted). 



to a given accuracy on a spatial lattice of sites. To this end we restore the lattice 
distance a, L = iVa, and write Eq{A)IL = a~^/7v(aA), So{A) = a~'^foo{aA), where 
£"0 does not include the counterterm. Figs. ^ and ^ show that LAEq = L{Eo{A) — 
EoiO)) = N\fNiaA)-h{0)) behaves like L\£oiA)-£om = N'{foo{aA)- f^iO)), 
modulated with 'oscillations' which in the scaling region turn into the physical LAEq 
(i.e. the quadratic (LA)'^ /it mod 27i potential). Since foo{aA) — /oo(0) = |r(aA)^ + 
T'{aAy + ■ ■ we expect that upon subtracting the counterterm, LAEq — ^t{LA)'^ 
equals the physical LAEq up to terms of order N'^{aA)'^ = N~'^{LA)'^. It follows 
that to given accuracy, LA can be at most of order ^/N. Hence the scaling region 
for LA grows as ^/N. 

In fig. 5 we display this scaling region for a lattice of size = 256. The 
nonperturbatively renormalized ground state energy LAEq — L'^ASq of the axial 
model (cf. fig. 4) is plotted superimposed on the gauge invariant (periodic) potential 
of the vector model. One can clearly see that several axial vacua near A = 
are degenerate and indistinguishable from the gauge invariant vector case, with 
violations of gauge invariance setting in at |C|^5. 

In the continuum limit iV — 00 we recover the equivalence with the vector 
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model. The scaling region then contains an arbitrarily large number of vacua |0)fc 



which transform into each other under large gauge transformations as in (f4.11|) , with 
the conserved fermion charge acting as 



Q|0), = -2A;|0)fc, 



(4.14) 



similar to (|4.12|) . Clearly, the conserved Q is not gauge invariant; it is the analogue 
of —Q5 in the vector model. 
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Figure 5: Scaling region for a lattice of N = 256 sites: the renormalized vacuum 
energy of axial QED is represented by the solid line, that of vector QED is dashed. 



5 Discussion 

We have shown using simple models in two dimensions, from spectral flow in the 
external field approximation that fermion creation and annihilation is possible even 
when fermion number is exactly conserved. The transition rate is determined by 
the usual anomaly formula, as if the fermion current had an anomalous divergence. 
We see no reason why the same resolution of the electroweak U(l) problem would 
not apply in four dimensions. 
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Furthermore, our study of the ground state energy of the axial model reveals that 
a standard perturbative mass counterterm for the gauge field, which was needed for 
restoration of gauge invariance, is insufficient nonperturbatively. Additional terms 
are needed ('oversubtraction') to get the global minimum of the ground state energy 
in the desired scaling region. In the vector model, the elegant way in which the 
staggered fermion method is able to give a bottom to the Dirac sea in a gauge 



invariant way (thereby circumventing the 'infinite hotel' explanation of the anomaly 
[0) is also worth noting. 

To complete our arguments, they have to be extended to the full theory with 
dynamical gauge fields. We also have to recover zero modes, fermion number viola- 
tion in correlation functions, etc.. Since we used continuous real time in the spectral 
flow arguments this may not be a problem. However, with discrete time e.g. in 
the euclidean formulation exact zero modes are non-generic and the consequences of 
this will have to be understood (for analogous questions in QCD see ref. [^, pU]). 
These questions lie outside the scope of our present inquiry. Nonetheless we shall 
briefly outline some expectations for the full theory with dynamical gauge fields. 

Consider first the vector theory as defined by the lattice regulated euclidean path 
integral in the compact formulation, in which the dynamical gauge field variables are 
phase factors t/^^r = exp(— zA^j,^) with A^j^^t ^ ("^i"; ^i"]; in lattice units. The path 
integral is then well defined without gauge fixing and its Hilbert space interpretation 
leads naturally to temporal gauge quantization, in which the integration over the 
time components U2xt of the gauge field provides the projection onto the gauge 
invariant subspace, the implementation of 'Gauss's law'. In the previous sections we 
used the Coulomb gauge however, so let us fix to this gauge: Uixt = Ur = exp(— zA,-) 
independent of x. In this case the integration over the time components U2xt in the 
path integral can be split into a summation over discrete large gauge transformations 
flkr- = exp{i2nkrx/L) plus fluctuations. The summation Pq = J2k^k over the large 
gauge transformations leads again to a projection onto the gauge invariant subspace 
whereas the integration over fluctuations leads to the Coulomb interaction. This 
interpretation of the integration over U2xt is appropriate in the scaling region. 

To compare with analysis in the continuum and also with the axial theory to be 
discussed below, let us assume that in the scaling region there is a which com- 
mutes with the hamiltonian (but not with Pq). This is not entirely obvious, since 
with fluctuating gauge flelds the spin-flavor interpretation of staggered fermions 
emerges only in the scaling region and even though A is constant, the local charge 
densities in the Coulomb interaction can only be expected to commute with such 
a Qs in the scaling region. In the scaling region the Hilbert space before pro- 
jection with Pq may then be split into sectors characterized by the eigenvalue of 
Qb: Qsl^)/; = 2fc|\E')fc, similar to the external gauge fleld case. Under large gauge 
transformations these sectors transform into one another, = and con- 
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sequently physical states Po|^)fc = Y^i have no well defined value of Q5. In 
particular Q5 does not annihilate the ground state |0), which is like a ^-vacuum 
superposition with 9 = 0, |0) = J2k |0)fc- 

The axial case is more difficult since we have not given a complete nonperturba- 
tive construction of the model with quantized gauge fields. The difficulty is the lack 
of gauge invariance at the fermion regulator scale. However, as mentioned in the 
Introduction, we have in mind an implementation which is gauge invariant on the 
scale of the lattice distance for the gauge field, which is much larger than the lattice 
distance for the fermions. Thus we may expect to regain invariance under gauge 
transformations which are smooth on the fermion regulator scale, in particular large 
gauge transformations, as we have seen in the scaling region of the external field 
model. 

This leads us to the conclusion that also in the axial case with dynamical gauge 
fields |0) = J2k |0)fc as well. Here however it is fermion number Q, which is anal- 
ogous to Q5 in the vector case, that is not well defined in the vacuum. Since the 
corresponding global U(l) symmetry is an exact symmetry of both the action and 
the fermion measure in the path integral, this suggests that the symmetry is broken 
spontaneously via the ^-vacuum: Q\0) = — I]fc2fc|0)fc 7^ 0. This picture though, 
needs to be reconciled with the expectation that spontaneous symmetry breaking 
occurs only in the infinite volume limit and that spontaneous breaking of a continu- 
ous symmetry in two dimensions is not possible. This will have to await the detailed 
construction of the model with dynamical gauge fields. 

Despite the simple nature of the models studied here, we hope the above picture 
is a valid description of the physics, also in similar four dimensional models where 
it mimics old ideas on the spontaneous breaking of U(1)a in continuum treatments 
of QCD. Our conserved vector current j'^ in the axial QED2 model is similar to the 
analogous axial current in QCD, = itpj^^'j^tp — 2n/C^ {uf is the number of fia- 
vors), which is conserved but not gauge invariant. The U(l) symmetry generated by 
the conserved charge Q5 = J d?x is in this scenario supposed to be spontaneously 
broken due to instanton — ^-vacuum effects. (There appears to be no universal agree- 
ment on this point p7| , pH]); the U(1)a problem in QCD can also be solved without 
invoking a gauge variant conserved jg |2^; for lattice implementations see [^.) 

To someone working exclusively in the continuum formulation, where models are 
constructed from formal unregulated expressions proceeding to well defined regu- 
lated expressions 'along the way', the gist of above remarks may appear so familiar 
that the question may arise, what's new? We repeat however, that with a non- 
perturbative regulator such as the lattice, models are well defined from the start 
and we have to determine the fate of the exact global U(l) symmetry and fermion 
number conservation. We have presented a resolution, at least in the external field 
approximation. 
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Appendices 

A Staggered fermion interpretation and ground 
state energy 



Various routes to the staggered fermion actions (|4.2| ) and ( [4.3|) from the continuum 
actions ( ^.6| , |2.2| ) have been described in refs. [Q, ^ jTU], Here we shall briefly 
review the continuum interpretation of the staggered fermion actions. 

The interpretation is done in momentum space on a lattice of N1N2 points; 
for simplicity we use lattice units a = 1. The momenta in the Brillouin zone are 
parameterized as + vr^^, with — 7r/2 < < +7r/2 and tca = {ttq = (0,0), tti = 
(7r,0),7ri2 = {n,n),7i2 = (0,7r)}, A = 1, ...,4. The regions of the Brillouin zone 
labeled by A provide the spin-flavor components of the fermion field according to 
il^p = Z]xexp[— i(p + ha)x]xx- For the simple case of constant A^, the actions 
(|4.2| , ^4.3|) can then be rewritten in the form ^ 

4 — _ 

for the vector case, and 

4 — / N 

5"/ = - Y.^P (^r^sinp^cosA^ + ^e^^r^cospi.sinA^)V'p, (A. 2) 

I\ll\2 



for the axial case. Here the tilde on J2p indicates the restricted momentum range 
given by = (n^ - l/2)(27r/A^), = -NJA -t- 1, • • • , A^/4. The 4 x 4 real 
symmetric matrices (r^)^^ come out ^ naturally as tensor products of Pauli 



matrices cr^, and r^. They commute with 4x4 real symmetric flavor matrices 
which correspond to shift symmetries of the staggered fermion action. The and 
are given by 

Fi = 0-3 , F2 = CTlTg , Si = (Tgri , = Tg . (A. 3) 

By a suitable unitary transformation the F's and S's can be brought into block 
diagonal form, F^,S^ l^l,^u, with 71,2,5 = cri,2,3, 6,2,5 = t-i,2,3- The classical 
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continuum limit corresponds to small p and A where the sines and cosines take 
their limiting forms. In the general case of non-constant gauge fields the classical 
continuum limit involves small momenta and values of A^^, as well. 

The coefficient r of the counterterm of the axial model can be obtained from the 
lattice vacuum polarization diagrams along the lines in ref. 0, which gives 



- + Ji + /2 ) ^ 0.3634 

TT 



2 



'^/^ d^q sin qi 1 



n/2 vr sin^ qi + sin^ q2 2 ' 
/, = i ^ (sin^g2-sin^gi)cosJg2 ^ _o^^g^7_ ^^ 4^ 

2 J-tt/2 TT L^^2 _^ gjj^2 



This is given in ref. for a somewhat different but equivalent action. In this 
sense the action studied in has less gauge symmetry breaking than the 'canonical 
action' ( [4. 3D studied here (notice that the value in ref. ^ had been normalized to 
one flavor). Below we shall obtain r again from the energy density of the ground 
state. 

The above representations ( [A.l| ) and ( |A.2| ) of the actions are useful for calculating 



the ground state energy Eq = MmN^^oo — {I/N2) Indet D, where D^y is the fermion 
matrix defined by Sf = —Xx^xyXy'^ we shall do so for the case A2 = 0, Ai = A 
considered in this paper. In momentum space det D = Opdet Dp, where Dp = 
iT2 sinp2 + ' ■ ■ can be read off from ( |A.1| ) and ( |A.2| ). We write Dp = T2{i sin p2 + TCp), 



where Tip is a hermitian matrix. Since det r2 = 1, the determinant of D is given by 
the product of the eigenvalues of {ismp2 + 'Hp), which come in complex conjugate 
pairs. Hence, for the evaluation of In det D we may use hi{zz*) = In^; + Inz*, with 
the principal value for the logarithm (In^ = ln|2;| +iaxgz, —it < argz < vr). For 
the vector theory we find 

Hp = T^sinip-A) (A.5) 

^ j'TT/2 dp2 

Eq = -Y. X!l^(^si^^'2 + xsin(p- A)), (A.6) 

p J-^/2 TT ^ 

and for the axial theory 

Tip = Fs cosAsinp + siny4cosp2 (A. 7) 

^ rTT/2 (ip2 ■ 

Eq = — / y^ln(isinp2 + XcosAsinp + siny4cosp2)5 (A. 8) 

where we have let A^2 — 00, p = pi, N = Ni, and the chirality x = =tl is the 
eigenvalue of F5. 
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As an example we describe the calculation for the axial case. The integrand 
has period vr, so we may replace jZl^/2^P'i/'^ ~^ /-7r^P2/2vr. The 'fermion doubling' 
introduced by this replacement is compensated by the factor 1/2. Next we change 
variables to z = exp ip2 and use 

Re/ -^ln(z-x) =ln(|a;|)^(|a;| -1), (A.9) 

J\z\=l ZTllZ 

where x is real and ^ is a Heaviside step function, 9{x) = 1, a; > 0, 9{x) = 0, x < 0. 
We factorize isinp2 + xcosAsinp + sinylcosp2 = (1 + sinA)(2; — z+){z — z_)/2z, 
with 

z± = ±exp(-e±^), e±x = ieo^ + In y^(l + sin A)/(l - sinA) , (A.IO) 
= ln(xs + Vs^ + 1) , sinheox = X'5 5 s = sgn(cos A) sinp. (A. 11) 

From the 'inverse Wick rotation' p2 = ip^ we see that z+ and e+ correspond to 
'normal fermion' excitations, whereas p2 = ip^ + vr shows that Z- and e_ correspond 
to 'doubler fermion' excitations introduced by doubling the momentum interval. 
Furthermore, the doublers have effectively opposite chirality, e„ = e+,-x- The step 
function in ( |A.9|) allows only the negative e±^ to contribute, as expected for the 
Dirac sea. The axial ground state energy can be evaluated as 

= A^ln2 - iVln(l + sin A) + 2^^ ^x^i'^x) > (A- 12) 

p X 

where = depends on A and p according to ( |A.11| ). 

The corresponding formulas for the vector case can be obtained from the above 
by letting A —>■ and then p p — A: 



= ln(x sin(p — y4) + y sin^(p — A) + 1) sinh e^^ = x sin(p — A) , (A.13) 

Eo = Arin2 + 2£^6^^^(-eJ. (A.14) 

p X 

Using e_(p, y4) = e+(p + vr,/!) this can be rewritten as a summation involving e+ 
only, over the whole periodic region p G (— vr, vr). This shows that Eq{A) is periodic 
in A with period 27r/iV. In the infinite volume limit the momentum summation may 
be replaced by an integration over p G (— vr, vr), and by p-translation invariance the 
energy density becomes independent of A, as expected from gauge invariance. 

Reintroducing the lattice distance a hy A ^ a A, Eq (iEq, L = Na, the 
continuum limit of the finite volume Eq can be found by expansion in a, taking A 
in the first period about A = 0. From e_(p, A) = e+{—p, —A) and the symmetry 
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of the momentum summation, only even powers of A occur. Using de-^-/dA 
— cos(p — A)/ coshe+ and replacing afterwards J2p by an integral for a 



Eo{A)-Eo{0) = LA' I ^ 

'-tt/2 it 



LA' 







dp d 



cos p smp 
+ sin^p)'^/^ 
cosp 



smp 



'1 + sin p) 



1/2 



+ ■ 



7r/2 TT 9p (1 + sin 



+ 



-LA' [1 + 0{a'A')] 



TT 



0, gives 

(A. 15) 
(A.16) 
(A.17) 



The two terms in (|A.15|) correspond to the 'ordinary' and 'leaf diagram' terms in the 
vacuum polarization diagram, and the differential identity in ( |A.16|) is essentially a 
Ward identity |11|. The factor two in comparison to the continuum formula ( |2.13|) 
is due to the two staggered flavors. 

Returning to the axial case, for which Eq is not gauge invariant without the 
counterterm, we write Eq = Nf]\f{A) (in lattice units). In the infinite volume limit 
the energy density Eq/N approaches a smooth function of A, 



fooiA) 



ln2-ln(l + sinA) +4 — e+{p, A)e{-e+{p, A)) 

J-tt/2 71 



-tt/2 TT 



(A.18) 



4 /-^ 

In 2 - ln(l + sinA) + - / 

TT J-n 



/4 COS 6 



' 1 + sin A 1 + sin 6^ 
1 — sin A 1 — sin 9 



where in the second line we used the substitution sinp = tan 6* and assumed A> 0. 
From foo{A) = \tA^ + 0(A^) we get an analytic expression for the coefficient of the 
counterterm, r = 1 — 2/7r = 0.363380... 



B Transfer operator and energy spectrum in vec- 
tor QED2 

In this appendix we construct the transfer operator for vector QED2 and rederive 
from it the energy spectrum. To the path integral a quantum mechanical Hilbert 
space is associated according to the coherent state formalism A rudimentary 

set of formulas for this purpose is given by 



\a) 
ai\a) 



exp(a|aj)|0) , {a\a) = exp(a^ai) , ^ = 'W j (ia^rfaj exp(— a^ai)|a)(a| 



aj|a) , {a\d\ = af {a\ ^ (a| exp(a|Mjjaj)|a) 



exp[aj {e 



ijCLj] 



(B.l) 
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In these appendices we distinguish operators in Hilbert space by a caret The di 
and a|, are a generic set of fermionic creation and annihilation operators with the 
usual anticommutation rules and 'empty state' |0), the and af are 'Grassmann 
variables', and M is some matrix. 

In the temporal gauge U2^x = 1, we can rewrite the vector action ( |4.2| ) in the 
form, using Xx = xtV2x, 

N2-1 

^ = - E Xr+l - Xt+lXr) + Xt^tXr] , (B.2) 

r=0 

where we have used (xi,X2) = {x,t) and suppressed the summations over x. The 
identification Xx,N2 = ~Xx,o, Xx N2 ~ ~Xxo takes into account the antiperiodic 
boundary conditions. The staggered fermion Dirac hamiltonian 7i is given by 

'Ht xy ~ 2 {jl2 x,T 1]l x,T x,T ^y,x+l + rj2 v,T rji y,T Uly,r ^x,y+l) , (B.3) 

and is manifestly hermitian. 

The time derivative terms —xtXr+i + xt+iXr are associated with the normaliza- 
tion factor exp(a"'"a) and measure factor exp(— a+a) in the coherent state formalism. 
This leads to the two time-slice construction described in refs. [|^, ^ |32|, where we 
interpret the x fields at odd time slices as 'type a' and at even time slices as 'type 
a"*"', writing 

Xr = \/2afc , xt = \/26fc , r = 2A; + 1, 

Xr = V2bt, xt = ^4, r = 2k, (B.4) 

where k = 0, . . . , N2/2 — 1 denotes pairs of time slices. Using these new variables 
the partition function can be rewritten in the form 



Z = Yl da^ dakdb^ dbk {ak+i,bk+i\f\ak,bk) (B.5) 

= Trf^2/^ (B.6) 

with 

{ak+i,bk+i\f\ak,bk) = e-^^^''^ e^+i'^fe+fc^+i^fe ^-2a+^,n2,+2b++, ^-2b,n2,+ia, _ ^g^^^ 

Recall that = A^^i is the total number of spatial sites. For the transfer operator T 
associated with the matrix element {ak+i,bk+i\T\ak,bk) we obtain the expression 
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The staggered transfer operator can be written as the square of the shift operator 
in the time direction p4|, but we do not need this here. 



We shall now determine the eigenvalues of T for the case Uix,t = exp(— iA) 
needed for the spectral flow {A may be considered constant in this calculation, as it 
is independent of the euclidean time r). First we consider the Dirac hamiltonian Ti. 
Similar to the two dimensional case in (|A.1|), the spatial Fourier transform brings 



7i in a more intelligible form 

E e-*(P+-^)^+^(''+-^)^ n.y = ia2)AB sin(p ~ A) N 5p,, , (B.9) 



= dAp, etc. (B.IO) 



X 



2 — - 

bHd = —^^^bpa2sm{p — A)dp, etc. (B.ll) 
p 

Here the indices A, B take only two values, tia = ttq, tti, and in a2 we recognize the 
upper block of = a^r-i- We will recover the complete Fs below. The eigenvalues 
of li, are given by A = ± sin(p — A), where ±1 is the eigenvalue of 02- The transfer 
operator reduces to Ilmodes 6xp(— 2 ln2) exp(— 2d"''A6^") exp(— 26Aa). 

For the moment we continue the calculation for a single mode. We seek linear 
combinations of a, 6^ and a^,6 which have simple commutation relations with T. 
This leads to a 2 x 2 real symmetric eigensystem for the coefficients of the linear 
combinations. From 

it is straightforward to check that 

n = (''J^ J,^,^i>f^e'-'^^^, (B.13) 
ln( A + VA2 + 1) , sinh 00' = X , (B. 14) 



with the orthogonal matrix V given by 

We can now retrace our steps and restore the matrix a2 to the eigenvalues of Ti in 
this basis: A — > (J2 sin(p — A), uo' ^ 0-2^^, turning ip^ into two component fields on 
which the a's act, hence obtaining four component fields ip on which both the cr's 
and r's act. Then the 2x2 orthogonal V gets replaced by a 4 x 4 unitary matrix 
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V, 

[hi) = (B.16) 
fiP = e^'^^'^^^f , sinhw = sin(p- A), (B.17) 

V = ^^^(e^^^^/^ + zrse""^/^) , (B.18) 
V 2 cosh to 

and we have recovered Fs = o"2T3- 

Defining the no-quantum state |0) by a|0) = 6^|0) = 0, it follows that ip\^) = 
and the eigenstates of the transfer operator are obtained by repeated application of 
ip'^ on |0), 

fV't|0) = |0), (B.19) 

The mode eigenvalues are obviously given by x^, with x the eigenvalue of T^. 
Under a large gauge transformation 

Q^a^Q = e*2-^/^a^ , , Q^blQ = e^^-^/^fet ^ (b.21) 
the ipp{A) transform as 



Q^p{A + 27r/N)Q = ^p_2VJv(^), p = Pmm + 27r/Ar, ■ ■ ■ ,p„ 
n^p^.JA + 2n/N)Q = T.E^iJp^^M), (B.22) 

where we used ai = T2'^2 and pmin = — Pmax = —7^/2 + tc/N. We see that the 
minimum momentum mode fiips its chirality x and fiavor ^ (eigenvalue of S5) as it 
changes into the maximum momentum mode, 

n^^^±{A + 2n/N)n = (-l)(x+0/2 ^-xr-^(A) . (B.23) 

The eigenvalue of FsSs remains the same. 

We shall now obtain operator expressions for the charges Q, Qe and Q5. Since 
fermion number Q = —i J2x J2 x,t is time independent we choose r = = even, and 
associate with Q = |(x^^2,TXr+i + Xr+i^^l rXr) = cl'^cl + bb^ the operator 

Q = a^a- Pb = a^a + W - N = ^Y.\^I^p ' ^ ■ (B-24) 

p 

This has the expected form (|3.3| ) for two fiavors. The other conserved U(l) charge 
Qe = -iJ2xj2x,t = Ex i(-exXrXr+i + ^xXt+iXr) IS also time independent (recall 
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ex = {—ly^ ^^). Choosing r = again leads to 

2 — 2 " — - 

iV p p 

= ^El^JrsSsi, (B.25) 

i V p 

with FsSs = —aiTi the flavor non-singlet axial generator corresponding to e^,. 

Although the chiralities x appear naturally in the mode spectrum of the stag- 
gered formulation, an axial current js^ or charge operator is somewhat extraneous 
to the formulation. A natural candidate is given by js^ = —ie^i,ju, Q5 = J2xjix,T- 
However, for a reasonable interpretation in staggered fermion theory, operators usu- 
ally have to involve some time average as well as space average. This leads us to 
investigate Q5 = | x,t + ji x,t+i) at r = 0. We find 

2 — 1 

<35 = -^I]2*^^P^2&p"^^^2a)cos(p- A), 

= -^E|V^;V^^ria2V^V^pCos(p-A), (B.26) 



2 



p 



r 



Tf'^Wp ] hTiSitanhLJ ^pCOs(p-A) 

A* „ ^ ^ Lcoshcij 



where FiSi = Ti, note that u depends on p — A according to eq. ( [B.17] ). In the 
scaling region cos(]9 — A) 1, u ^ and the rather ugly looking Q5 reduces to the 
continuum form. Away from the scaling region there is a reduction of the 'strength 
of Fs' in Q5 and flavor symmetry breaking. 

The operators Q, and Q5 are invariant under the large gauge transformations 
( [B.22| ), as is obvious from their original gauge invariant expression in terms of the 
staggered fermion fields. We can define similarly a conserved chiral charge operator 

Q5 which commutes with the transfer operator ( [B.8| ) but is not gauge invariant 
under the large gauge transformations (|B.22D , 



^5 = ^ElV^Jrs^p. (B.27) 



The chiralities x = =tl of the modes ipR^L = |(1 i Ts)^ are the eigenvalues of this 
operator Q^. 
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C Transfer operator and energy spectrum in axial 
QED2 

To calculate the transfer operator and its eigenvalue spectrum in the axial-vector 
model we follow the same steps as in the previous section, and will therefore only 
point out the differences in the calculation as compared to the vector case. Using 
the definitions ( p.4|) , we find for the action (|4.3| ) in terms of the fields Xr and xty 
the expression 

^-2-1 _ 

S = -T.[U^XtXr+l~PxUlXr)+xtnrXr], (C.l) 
r=0 

a = 1 + sinA, /3 = l-smA. (C.2) 

The matrix Ti-r xy in is given by the same expression as in ( [B.3| ), except that 

Ui x,T has to be replaced by the real form cos A. With the same assignment of the 
variables ak,-.-,b^ to the Xt and xt fields, as in ( |B.4[ ), and after the rescaling 
Ofc (Ofc) — >■ ttk/y/a (a^/i/a) and bk (fc^) — > bk/y/P (b'l^/VP) we find that the matrix 
element of the transfer operator is given by 

{ak+i,Ok+i\i \ak,Ok) — e e 

xe-24+i(^2fc+2/v/a^)''fc+i g(/3/"X+iafc+W/3)fe;t+ifefc g-2fefc{W2fc+i/v^)"fc . (C.3) 

The second exponential factor in (|C.3| ) arises from the jacobian in the above rescaling 
of the Grassmann measure: d{ak/\fa) = ^Jadak, etc. Using the rules ( |B.1| ) for 
transcribing the Grassmann algebra matrix elements to their associated Fock space 
operators we obtain from ( |G.3| ) 

^ g-2Afln2 g+Afln(a/3) ^-2~a'^Hoh'^ gln(/3/a) (ata-fotfo) g-2feHoa (^Q 4^ 

Horxy = Sgn (cOsA)i (?72 x,t ^71 a;,T ^^^.ai+l + ^2 y,r ^71 y.r • (C.5) 

For sgn cos A > 0, which is satisfied because we require \A\ < 7r/4 (cf. sect. 4), the 
matrix TYq is just the free staggered fermion Dirac operator. In the factor d^d — wb 
we recognize the charge operator Q of ( |B.24| ). The factor exp[ln{a / f])Q] has simple 



commutation relations with the fields d,¥ and d\b and the analogue of ( p.l8|) in 
this axial case is given by 

fijj = eMa//3)+2<^o^2r3 ^ sluh c^q = slu p , (C.6) 

V = , ^ (e-^o<^2/2 ^ ^^^e'^o'X2/2^ _ 
V2 cosh tuo 

The interpretation of the operators -j/^, follows that of the vector case. The energy 
spectrum following from ( |C.6| ) is the same as in ( [A.10| ). 
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